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Cosmology and short baseline neutrino oscillation data both hint at the existence of light sterile 
neutrinos with masses in the 1 eV range. Here we perform a detailed analysis of the sterile neutrino 
scenario using both cosmological and SBL data. We have additionally considered the possibility 
that the extra neutrino degrees of freedom are not fully thermalised in the early universe. Even 
when analyzing only cosmological data we find a preference for the existence of massive sterile 
neutrinos in both (3+1) and (3+2) scenarios, and with the inclusion of SBL data the evidence is 
formally at the 3.3a level in the case of a (3+1) model. Interestingly, cosmological and SBL data 
both point to the same mass scale of approximately 1 eV. In the (3+1) framework WMAP9+SPT 
provide a value of the sterile mass eigenstate 7714 — (1.72 ± 0.65) eV: this result is strenghtened by 
adding the prior from SBL posterior to m4 = (1.27 ± 0.12) eV (m4 = (1.23 ± 0.13) eV when SDSS 
is also considered in the cosmological analysis). In the (3+2) scheme, two additional, non-fuUy 
thermalized, neutrinos are compatible with the whole set of cosmological and SBL data, leading to 
mass values of m4 = (0.95 ± 0.30) eV and = (1.59 ± 0.49) eV. 



I. INTRODUCTION 

The past decade has seen a rapid increase in our un- 
derstanding of neutrinos. Oscillation experiments have 
now established the main structure of the leptonic mixing 
matrix and provided evidence for at least two neutrino 
mass states of non-zero mass (see e.g. Ref. [T] for a recent 
overview). However, some crucial questions are still left 
unanswered: the absolute mass scale of neutrinos is ex- 
tremely hard to measure in laboratory experiments and 
is therefore currently only poorly known; the hierarchy 
of neutrino masses is not yet disentangled; the possible 
existence of additional neutrino species beyond the three 
predicted by the Standard Model is still an open issue. 
In fact, short-baseline (SBL) oscillation experiments, as 
well as reactor neutrino flux measurements seem to hint 
at the existence of a fourth (3+1 models, possibly with 
non standard interactions) or fifth (3+2 models) sterile 
neutrinos, mixing with the active neutrinos and with a 
mass around 1 eV, see Refs. [2H2] . 

Cosmology has, at the same time, provided impor- 
tant insights into some of these questions. Neutrinos are 
produced in copious amounts in the early Universe, and 
are still present as a cosmic neutrino background. Even 
though this background is extremely difficult if not im- 
possible to measure directly, it influences a number of 
cosmological observables such as the cosmic microwave 
background (CMB), and the power spectrum of matter 
fluctuations. The effective number of relativistic degrees 
of freedom changes both the position and the shape of 
peaks of the CMB temperature power spectrum at high 
multipoles "W . The neutrino mass affects both the CMB, 
through the enhancement of the early integrated Sachs- 
Wolfe effect, and the matter power spectrum, via the 
free-streaming that suppresses the power on small scales. 
Thanks to these fingerprints, cosmology can strongly con- 
strain the absolute neutrino mass scale and the cosmic 
neutrino background [SI HUj . 



In the past years cosmology has provided some hints for 
a non standard value of the effective number of relativis- 
tic degrees of freedom, see Refs. [TTlI14j . with the pre- 
ferred value in the late-time universe (around, or subse- 
quent to, recombination) being significantly higher than 
the standard model prediction of A'^eff = 3.046 [U]. Such 
additional relativistic energy density is usually referred 
to as dark radiation. The light sterile neutrinos hinted 
at by SBL data would be an excellent candidate for dark 
radiation [TB], even though misinterpretations of the na- 
ture of this non standard A"eff can arise from degeneracies 
between N^s and other cosmological parameters [IT. 

Recently, new CMB data released by the WMAP col- 
laboration [18], the South Pole Telescope (SET) [19], and 
the Atacama Cosmology Telescope (ACT) [20] have led 
to a somewhat confusing situation concerning dark ra- 
diation and, in general, those neutrino properties that 
can be constrained by cosmology: the WMAP9 and SPT 
data both confirm the presence of an extra dark radia- 
tion component, while the new ACT data seem pointing 
towards a value of A^cfi in agreement with the standard 
model prediction, in contrast with their own previous 
analysis [Hj. 

Fortunately this problem might be resolved by new 
data from the Planck mission [22] which should be able 
to constrain Ngg with a much better precision than ex- 
isting data (see e.g. Ref. [531 HI])- However, even if 
the incoming cosmological data from Planck confirms the 
standard cosmological value for A^cff, cosmology should 
address the discrepancy with the SBL neutrino oscilla- 
tion results. Since the SBL experiments do not provide 
direct constraint on the cosmological thermalization of 
the sterile neutrinos, a small lepton asymmetry can be 
invoked to reduce the thermalization efficiency (see e.g. 
Ref. [^ [55] for recent treatments) and reconcile these 
two branches of neutrino physics. 

So it is indeed timely to investigate the interplay be- 
tween neutrino oscillation experiments and cosmology in 
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determining the neutrino properties [571 HH] • This paper 
is aimed at investigating the joint constraints on neu- 
trino number and neutrino mass from these two different 
branches. The paper is organized as follows: in Section 
[n]we present the SBL analysis data, which includes up- 
dates on our previous analysis; in Section [TTT| we describe 
the cosmological data used in this work and the method 
we have applied to analyze them; in Section |IV| we pro- 
vide an update of the cosmological constraints on sterile 
neutrinos; in Section |V] we perform the joint analysis by 
applying a short-baseline prior to the cosmological anal- 
ysis. Finally, Section \VJ\ contains our conclusions. 

II. SBL ANALYSIS 

We consider 3+1 and 3+2 extensions of the standard 
three-neutrino mixing (see Ref. [29^) in which we have 
one (7714) or two (7714 and 7715) new neutrino masses at the 
eV scale and the masses of the three standard neutrinos 
are much smaller: 

mi, m2, TO3 < m4 < ?7i5 . (1) 

In this case, the squared-mass differences 

Amli ~ ml and Am^^ ~ m^ , (2) 

with Am^j = ml — to|, generate short-baseline oscilla- 
tions through the mixing relation 

4 or 5 
fc=l 

between the flavor neutrino fields Va (ct = e, /i, r, si, S2 in 
the 3+2 model) and the neutrino fields with masses 
mk (fc = 1, . . . , 5 in the 3+2 model). U is the unitary 
mixing matrix. 

In the analysis of short-baseline data we adopted the 
approach described in Ref. 127] , with some improvements 
in the considered data sets, which can be divided in the 
following three groups: 

1. The — >■ Ve and — >■ appearance data of the 
LSND [3^, KARMEN NOMAD [32], Mini- 
BooNE [33] and ICARUS [SJ experiments. In par- 
ticular, we use only the MiniBooNE data above 475 
MeV, because the data at lower energy contains an 
anomaly which cannot be explained by neutrino os- 
cillations [HH] (an interesting possiblity of reconcil- 
ing the low-energy anomalous data with neutrino 
oscillations through energy reconstruction effects 
has been suggested and discussed in Refs. [351135] ). 

2. The z/g and De disappearance data described in Ref. 
[37] . which take into account the reactor [5S] and 
Gallium ^39] anomalies. 

3. The constraints on i/^ and disappearance ob- 
tained from the data of the CDHSW experiment 





3+0 


3+1 


3+2 


Xmin 


280.2 


236.1 


229.0 


NDF 


230 


227 


223 


GoF 


1.3% 


32% 


38% 


Ami, [eV^] 
|(7e4p 




1.62 


0.91 




0.031 


0.015 


Ami, [eV^] 




0.01 


0.011 
1.61 
0.0226 








0.00664 


V 






1.567r 


Axpg 
NDFpG 




6.6 

2 


11.12 

4 


GoFpG 




4% 


2.5% 



TABLE I. Values of Xmim number of degrees of freedom 
(NDF), goodness-of-fit (GoF) and best-fit values of the mix- 
ing parameters obtained in our 3+0, 3+1 and 3+2 fits of 
short-baseline oscillation data. The last three lines give the 
appearance-disappearance parameter goodness-of-fit (PG). 

[40] . from the analysis [41] of the data of atmo- 
spheric neutrino oscillation experiments, from the 
analysis [3] of the MINOS neutral-current data 
[32] and from a new analysis of the SciBooNE- 
MiniBooNE neutrino j43, and antineutrino [44] 
data. 

The results of a least -squares analysis of the SBL data 
is presented in Tab. [I] First, we notice that the values of 
the goodness-of-fit obtained in the 3+1 and 3+2 models 
are satisfactory and much better than what is obtained in 
the case of absence of sterile neutrinos (3+0). We notice 
also that the value of the appearance-disappearance pa- 
rameter goodness-of-fit is acceptable and it is remarkable 
that it is better in the 3+1 model than in the 3+2 model, 
contrary to previous results with older data [H [31 US] . 
The reason is that the values of Xmin ^'^^ appearance 
(APP) and disappearance (DIS) data in the two mod- 
els are: 

(xL„)l'tp = 50.4, ixLJms = 179.1 , (4) 

(xii„)it.p - 40.1, {xLnfuil = 177.8 . (5) 

Since the 3+2 model can fit significantly better than the 
3+1 model only the appearance data, the larger differ- 
ence between the global xii„ and (Xmin)APP + (Xmin)Dis 
in the 3+2 model cannot be compensated by the increase 
of the number of degrees of freedom. 

Figure p] shows the allowed 3cr regions (99.73% CL) 
in the sin 2'de^i~Am1i plane in the 3+1 model obtained 
with disappearance data, appearance data, and all data 
respectively. The quantity sin^ 2i?e^ = 4|[/e4p|C^M4l^ 
the amplitude of and — ?> v^, oscillations 

(see Ref. [2^). One can see that the combination of 
the constraints obtained from disappearance data ex- 
cudes the large-sin^ 2??e^ part of the regions allowed by 
appearance data. From the global fit we obtained a 
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FIG. 1. Allowed 3cr regions (99.73% CL) in the sin^ 21?^^- 
Am|i plane in the 3+1 model obtained from Ve and dis- 
appearance data (left of the dark red curve), and dis- 
appearance data (left of the dark green curve) , combined dis- 
appearance data (left of the dark orange curve), — >■ Ve and 
S'fi — >■ i'e appearance data (inside the blue curve) and from 
the global fit (inside the red curves). The best-fit points in 
the last two cases are indicated by crosses. 



low-Am|]^ allowed region at Am|]^ ~ 0.8 — 2 eV^ and 
sin^ 2'detj. — (0.5 — 3) x lO^'', containing the best-fit point, 
and a high-Aml]^ allowed region at Amli ~ 6eV^ and 
sin^ 2i?e^ ~ (0.8-2) x lO'^. 

The appearance-disappearance tension discussed in 
previous papers (e.g. Refs. [5H51 US]) is still present, 
because, as one can see from Fig. [l] the best-fit point 
of the appearance data is excluded at about Scr by the 
disappearance data. However, the tension is less severe 
than that which was obtained with old data, as testified 
by the acceptable parameter goodness-of-fit in Tab. |l] 

In the following we combine the results of the analysis 
of short-baseline and cosmological data considering both 
3+1 and 3+2 schemes, in spite of the fact that the 3+1 
scheme is sufficient to explain the current short-baseline 
data and the 3+2 scheme is disfavored by Occam's razor. 
We consider also the 3+2 scheme because of the current 
interest in it (see, for example, Ref. [46]) and because 
future data may reverse the preference. 



III. COSMOLOGICAL METHOD 

We have modified the Monte Carlo Markov Chain 
(MCMC) pubfic package CosmoMC gT] (October 2012 ver- 
sion) to account for the datasets listed below and in order 



to sample the extended parameter space for our (3+1) 
and (3+2) investigations. 

The cosmological analysis is performed by employing 
various combinations of datasets. The WMAP 9-years 
data release 18 represents our basic CMB dataset. At 
high multipoles, we use new data from the CMB experi- 
ments South Pole Telescope and Atacama Cosmology 
Telescope The information on dark matter cluster- 
ing comes from the matter power spectrum derived from 
the Sloan Digital Sky Survey Data Release 7 luminous 
red galaxy sample [15]. We have also investigated the 
impact on our constraints of the Baryonic Acoustic Os- 
cillations (BAO) results of Refs. [15HST] and of the prior 
on the Hubble constant coming from the Hubble Space 
Telescope measurements |52j . 

Short-baseline results have been included as a prior 
in the analysis. Therefore the final of the combined 
analysis is simply given by 



2 _ 2 
Xtot Xcosmology 



XSBL- 



(6) 



Concerning the inclusion of the SBL information in the 
cosmological analyses, the difference with respect to our 
previous work [37], is that here we directly incorporate 
the coming from SBL data (analyzed in a variety of 
frameworks) in the MCMC sampling: this allows for un- 
biased constraints on those parameters that are common 
to the cosmological and SBL analyses. 

Our basic cosmological model is the 6-parameters flat 
ACDM model. Its parameters are: the physical baryon 
and cold dark matter densities flijh^ and ri^j^/i^, the ra- 
tio of the sound horizon to the angular diameter distance 
at decoupling 9s, the optical depth to reionization r, the 
scalar spectral index rig and the overall normalization 
of the spectrum Ag. In the cosmological analyses, we 
also account for the Sunyaev-Zeldovich (SZ) effect and 
foreground contributions by including three extra ampli- 
tudes: the SZ amplitude asz, the amplitude of clustered 
point sources Oc and the amplitude of Poisson distributed 
point-sources Op. Since the ACT team uses a different 
mask for identifying and removing point sources, when 
the ACT data are included in our analysis (Section IV) 



a fourth extra amplitude is needed to account for the 
Poisson contribution in the ACT data apACT different 
from apSPT- Furthermore when the ACT likelihood is 
included in our analysis, we split the SZ contribution in 
two amplitudes referred for the two different SZ effects: 
o-ksz for the kinetic SZ effect and atsz for the thermal 
SZ effect. 

In order to include extra sterile neutrinos, the ACDM 
model is extended to a A Mixed Dark Matter model 
(AMDM) by introducing a hot dark matter component 
in the form of massive neutrinos. Usually this compo- 
nent is parameterized as the neutrino mass fraction 
(the neutrino-hot dark matter density over the total dark 
matter density), as in Section IV but in the joint analy- 
sis of Section |V] the masses of the single mass eigenstate 
are used as free parameters, in order to directly sample 
the parameter which the SBL prior acts on. 
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With the above assumptions, the contribution of mas- 
sive neutrinos to the energy budget of the Universe fol- 
lows the usual relation: 



93.14 eV 



(7) 



where ^ rrii, denotes an effective sum of masses, specified 
below for the different cases under study. 

Before attempting the combined (cosmological + SBL) 
analysis, we perform a number of tests on the cosmolog- 
ical data alone, in order to assess the requirements on 
dark radiation from the various cosmological data sets. 
These results will then be compared to the interpretation 
of dark radiation in terms of sterile neutrinos in the (3-1-1) 
and (3-)-2) models. From the point of view of cosmology 
there is no difference between active neutrinos and ster- 
ile neutrinos provided they are fully thermalized. For the 
typical mass differences and mixing angles hinted at by 
the SBL data this is indeed the case (see e.g. \^5j2Qj). In 
this case only the total sum of neutrino masses and the 
total (effective) number of neutrino species are relevant 
parameters (with the current level of precision of cosmo- 
logical data). However, it is entirely possible that the 
sterile states are only partially thermalized. This can for 
example happen in models with non-zero lepton asymme- 
try [251 HB]. Therefore, in Section IV we perform three 
different analyses: 

(A) All neutrinos are massless: in this case cosmological 
observables are just sensitive to the total effective 
number of relativistic (fermionic) degrees of free- 
dom denoted by Ncs, where A'^eff = 3.046 if 
only (effectively massless) standard model neutri- 
nos are present; 

(B) Neutrinos are allowed to be massive, with a com- 
mon total mass ^ to^ and a number effectively 
equal to N^ff- This corresponds to: 



A^: 



:ff X niM 



(8) 



(C) The three standard model neutrinos are (effec- 
tively) massless, while extra sterile neutrinos are 
present in (effective) number equal to Ns and total 
(effective) mass equal to ms- This corresponds to: 



rut^s = Ns xms 



(9) 



and to a total effective number of neutrinos Aoff — 
3.046 -I- Ng. This case differs from case (B) since 
here the presence of 3 massless neutrinos fully con- 
tributes to radiation energy, while in case (B) all 
neutrinos are effectively massive and there may 
be a different impact on matter-radiation equal- 
ity, CMB anisotropics and cosmic structures for- 
mation. Clearly in this case Ns measures only the 
extra-number of degree of freedom responsible for 
dark radiation, while in case (B) A^^g is the total 
number of (potentially) massive neutrinos. 



We have considered the possibility that sterile neutri- 
nos are only partially thermalised. In our analysis we 
express this in terms of the following parameter: 



A,; 



rij 



(10) 



where rii denotes the actual number density of a sterile 
neutrino of mass while nf^ denotes the number density 
of a standard neutrino (fully thermalized, with a Fermi- 
Dirac phase-space distribution) with the same mass m^. 
The multiplicity parameter A^^ therefore defines the frac- 
tional contribution of the (non-standardly thermalized) 
sterile neutrino to the dark matter energy density, and 
is a number defined in the interval [0, 1]. In our assump- 
tions, Ni is considered to be independent from the mass 
TO,;, and therefore A^^ and can be treated as indepen- 
dent parameters in the analysis. While this may not be 
the most general case, it is generic enough to study the 
possibility to allow, in the cosmological data, for extra 
sterile neutrinos endowed with the properties dictated 
by the SBL studies. Examples of values of A'i < 1 may 
be related to partial thermalization (possibly related to 
the presence of a lepton asymmetry, since the mixing 
angles obtained from the SBL analysis are large enough 
to ensure thermalization of the sterile neutrinos if the 
asymmetry is absent) or non-standard phase-space dis- 
tributions for sterile neutrinos. The limit Ni = 1 refers 
to a fully thermalized Fermi-Dirac sterile neutrino. 

In out joint (cosmological + SBL) analysis (discussed 
in Section |v]) we assume the three active neutrinos as 
massless and the J2 Eq. ([t]) is therefore given by 

^ = NiUii, where the index i runs over the num- 
ber of sterile neutrinos considered in our analysis: 

^m,, = A^4XTO4 (3+1) (11) 

^m,, = Ni X 7714 + N5 X "^5 (3+2) (12) 

We have considered a top-hat prior in the range [0, 1] 
for the multiplicity parameters, while the sterile neutrino 
masses have been subjected to a prior given by the 
results of the SBL analysis. 

Finally, in order to grasp a connection between the 
joint (cosmological + SBL) analysis and the more typ- 
ical cosmological investigations (where the relevant pa- 
rameters are just the sum of the neutrino masses and an 
effective number of additional light states, as discussed 
above), we have re-analyzed the (3+2) case by using the 
SBL-prior on 7714 and TO5 projected over the sum of the 
two masses TO4 + 777,5. This case corresponds again to a 
situation where: 



A"4 X 7774 + A^s X (3+2) (13) 



where we look at the results in terms of (A4 + A's) and 
7774 + 7775 . This case is discussed in Section |V] and differs 
from the case studied in connection with case (C) because 
information coming from the SBL studies is included in 
the analysis. 
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FIG. 2. Case (B): N'^^ massive neutrinos with total mass 
^ m„ - Two dimensional marginalized 68% and 95% confi- 
dence level regions in the plane J]] 'm^ vs N'^^ . Leftmost (red) 
contours refer to CMB-only data (WMAPO+SPT+ACT), 
while rightmost (blue) contours include also BAO and HST. 
The larger (green) area denotes the results for WMAP9+SPT 
datasets. 



IV. COSMOLOGICAL RESULTS 

The cosmological analyses performed on the full CMB 
dataset (WMAP9, SPT and ACT combined) is shown 
in Table llll The results for case (B) are also illustrated 
in Figure]2] where the marginalized 68% and 95% con- 
tours for N'^^ and ^ are reported. The analysis for 
the combination of WMAP9-t-SPT-hACT datasets refers 
to the leftmost (red) contours; the further inclusion of 
BAO-I-HST moves the contours to the right (blue con- 
tours), referring to larger iV^fj, accompanied by smaller 
values of ^TOj^). Figures shows the one dimensional 
marginalized posteriors for A'^g and ^m^. 

The results shown in Table and Figures|2]and[3]point 
toward the conclusion that the full set of present CMB 
data do not show any indication for a non-standard value 
of the effective number of relativistic degrees of freedom. 
Nevertheless this result is mostly driven by the new ACT 
likelihood and a strong tension actually emerges between 
the new SPT and the new ACT data. This tension turns 
out in a bias that has a strong impact on the cosmolog- 
ical results, especially on the neutrino sector. Removing 
the ACT data and using only WMAP9 and SPT, in the 
massless neutrino case (case (A)) we recover a 1.5cr evi- 
dence for an extra dark radiation component: 

AT^ff = 3.87 ± 0.55 (Icr) case (A) 

WMAP+SPT 

This result is fully consistent with those of the WMAP9 
[TB] and SPT [19 analyses. In case of massive neutrinos 
(case (B)) Figure [2] shows that, when all CMB data are 



considered (red contours), an anti-correlation emerges 
between N'^^ and the total mass ^ m^: a higher value of 
the sum of the masses seems to be consistent with a lower 
value of the number of neutrinos that share the same 
mass. This inconsistency can be traced to the tension 
between the high multipole datasets (ACT and SPT). In- 
deed, excluding the ACT data, we get the green conturs 
of Figure [2] the higher mass values allowed by SPT cor- 
respond to higher numbers of neutrino species, and there 
isn't any strong correlation between these two quantities. 
Furthermore the above tension between ACT and SPT 
is less pronounced also if we take into account BAO and 
HST data. In this case (blue contours) we recover the ex- 
pected positive correlation between 'Yli'^v and A^^ff, be- 
cause the HST measurements fix the expansion rate and, 
as a consequence, the dimension of the sound horizon at 
recombination and the damping scale. So an enhance- 
ment in the sum of the masses is needed in order to get 
a higher value of the number of neutrino species. More- 
over, considering WMAP9-FSPT-hACT-FBA0-f HST, we 
get a slight preference for a non-standard value of the 
effective number of relativistic degrees of freedom (see 
Figure |3|: 

A^^ff = 3.60 ± 0.35 (Icr) case (B) 

WMAP+SPT-fACT+BAO+HST 

Concerning the mass, the constraint on ^ is signifi- 
cantly tightened by adding the low redshift observations 
(BAO and HST) data and it turns out to be: 

^ < 0.59 eV (95%c.l.) case (B) 

WMAP+SPT+ACT+BAO+HST 

with a best fit value of ^ m^j — 0.23 eV, as we can see 
in Figure [3j 

Excluding ACT data and using only WMAP9-I-SPT, 
in case of massive neutrinos (case(B)), we obtain: 

N'^s = 3.44 ± 0.56 (Icr) case (B) 

WMAP-fSPT 

with a total mass 

^m^< 1.85 eV (95%c.l.) case (B) 

WMAP+SPT 

In Figure [s] the dot-dashed (green) line shows the one 
dimensional marginalized posteriors for ^ and A'^g 
in this case with only WMAP9-FSPT. The best fit value 
for the sum of the masses is shifted beyond 1 eV and the 
range of possible mass values is much less constrained 
than in the above cases when ACT data are included. 



V. JOINT ANALYSIS 

Since the new SPT [1^ and ACT 20J data seem to 
point towards opposite directions concerning the effec- 
tive number of relativistic degrees of freedom, here we 
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Parameters 




Case (A) 


Case (B) 


Case (C) 






WMAP9+SPT+ACT 


WMAP9±SPT±ACT 


WMAP9±SPT±ACT 






0.02259 ± 0.00045 


0.02224 ± 0.00053 


0.02281 ± 0.00035 






0.115 ±0.008 


0.120 ±0.009 


0.126 ±0.009 


Os 




1.0418 ±0.0013 


1.0417 ±0.0013 


1.0416 ± 0.0011 


T 




0.086 ±0.014 


0.086 ±0.013 


0.088 ±0.013 


Us 




0.972 ±0.017 


0.959 ± 0.020 


0.966 ± 0.010 






3.178 ±0.041 


3.213 ±0.049 


3.212 ±0.039 




(A) 


3.17 ±0.47 


— 


— 




(B) 


— 


2.97 ±0.48 


- 


Em.[eV] 


(B) 


- 


< 1.17 


- 


Ns 


(C) 


- 


- 


< 0.89 


Sm,s[eV] 


(C) 


- 


- 


< 2.39 


//o[km/s/Mpc] 


71.7 ±3.4 


64.9 ±5.5 


68.9 ±2.0 


0-8 




0.822 ±0.029 


0.714 ±0.073 


0.699 ±0.071 






0.269 ±0.020 


0.345 ± 0.063 


0.315 ±0.034 


dtSZ 




4.4 ± 0.9 


4.4 ±0.9 


4.5 ±0.9 


aksz 




< O.Z 


< 0.4 


< 6.1 


ttc 




6.0 ±0.5 


6.0 ±0.5 


6.2 ±0.5 


dpSPT 




18.5 ± 1.6 


18.6 ± 1.6 


18.4 ± 1.6 


UpACT 




7.0 ±0.3 


6.9 ±0.3 


6.9 ±0.3 


Xmin 


8962.2 


8962.9 


8961.3 



TABLE II. Values of the cosmological parameters and their 68% confidence level intervals for the the three cosmological analysis 
described in Section III case (A) refers to A^eff massless neutrinos; case (B) refers to N'^ff massive neutrinos with total mass 
'^rriv; case (C) refers to 3 massless (active) neutrinos plus Ns additional (sterile) neutrinos, with total mass mg. Upper 
bounds are quoted at 95% C.L. 




FIG. 3. Case (B): A^^g massive neutrinos with total mass '^rriy - One dimensional marginalized posteriors for N'^g (left panel) 
and 5i;m^ (right panel). Solid (red) hues refer to CMB-only data (WMAP9±SPT±ACT), while dashed (blue) hues include 
also BAG and HST. The dot-dashed (green) curve denotes the results for WMAP9±SPT datasets. 



adopt the conservative approach to do not combine the 
two datasets in the joint analysis with SBL experiments. 

Moreover, we decided to consider only the SPT dataset 
for the following reasons. First of all, the ACT dataset, 
when only the temperature angular power spectrum data 
is considered, provides an indication for gravitational 
lensing larger by 70% than what expected in the stan- 



dard ACDM model (with massless neutrinos) at more 
than 95% c.l. (see [2Qj but also [53]). Since massive neu- 
trinos decrease the gravitational lensing signal, the lens- 
ing anomaly biases the ACT results on neutrino masses 
towards more stringent constraints (see Figure 3, right 
panel of |53j). The origin of this anomalous signal is un- 
clear and could possibly be due to a systematic in the 
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data. The SPT dataset, on the other hand, exhibits no 
anomaly in the lensing signal. 

Secondly, ACT is composed by two maps in two differ- 
ent regions of the sky, defined respectively as the ACT- 
E and the ACT-S datasets. The ACT-S map overlaps 
with the region sampled by SPT and the two should not 
be combined because of their covariance. The ACT-E 
dataset (that constitutes approximately 50% of the en- 
tire dataset) provides constraints that are even in larger 
tension with SPT (see [20]). 

In what follows we report the results for the joint (cos- 
mological + SBL) analysis, obtained with different com- 
binations of cosmological datasets. The SBL posterior 
probabilities on the neutrino parameters are used as pri- 
ors in the MCMC sampling of those parameters relevant 
to cosmology, i.e. the sterile neutrino masses. In the 
(3-1-1) scheme we have a prior on the single extra neu- 
trino mass m4; in the (3-1-2) scheme we have priors for 
the two extra neutrino masses 7714 and m^; finally, in 
(3-1-2) scheme we also consider the only relevant parame- 
ter for cosmology, which is the sum of the neutrino masses 
m4 -|- m5 and for that we use a prior obtained from the 
SBL analysis by projecting the SBL priors for the single 
masses to the sum of the two. 



A. (3+1) Joint Analysis 

The results of the joint (cosmological -I- SB L) a nalysis 
for the (3+1) scheme are reported in Table III and in 
Figure |4] and [Sj The analysis refers to the case of Eq. 
(111. In this analysis, in addition to the CMB datasets 



sbl 

5.0" cosmology 

cosmology+sbl 



(WMAP9 [I8j and SPT pJJ) we also consider information 
on the matter power spectrum coming from large scale 
observables: we include data from the Sloan Digital Sky 
Survey Data Release 7 (SDSS) 

When only cosmological data is used, bounds on N4 
and 7714 are broad: the allowance for an extra neutrino 
which cosmologically acts only as a fraction of a fully 
thermalized one is compatible with cosmological data for 
a neutrino mass as large as 2.09 eV. The bound on the 
multiplicity parameter A^4 is as large as 0.96. These re- 
sults stand when also SDSS data are included: CMB- 
only data, instead, have a preference for lower values of 
A'4 and 771,4, with best-fit values different from zero at 
about the 3a level. The central value for N4 is 0.65 and 
deviates from zero at 2.9ct; the preferred value for m4 is 
1.72 eV, and differs from zero at the 2.7cr level. 

When the SBL posterior probability on 7714 in used as 
a prior in the cosmological analysis, the situations moves 
toward a clear preference for a neutrino mass above 1 eV 
and slightly lower values for the multiplicity parameter. 
The result on the mass is expected, since SBL data have 
a clear preference for a non-zero neutrino mass, and this 
reflects on the posterior probabilities of the joint analy- 
sis. The combination of CMB-only data with SBL results 
does not change appreciably the multiplicity parameter 
(preferred value of 0.69, different from zero at the 3cr 



^ 3.0- 




[eV] 



FIG. 4. (3+1) analysis - One dimensional marginalized pos- 
terior for 7714. The thick (green) and thin (blue) lines refer 
to the case of Table |III| with and without the SBL prior, re- 
spectively. Solid lines stand for the analysis on CMB-only 
(WMAP9+SPT) data; dot-dashed lines refer to the inclusion 
of information from the matter power spectrum. The (red) 
dashed line shows the 3+1 SBL posterior. 95% C.L. upper 
bounds on the mass for the different cases are reported as 
vertical lines. 



level) while it slightly reduces the best-fit value for the 
mass: 1.27 eV, with an error of 0.12. The inclusion of 
SDSS data enlarges the allowed interval for the multiplic- 
ity parameter, while the preferred value for the masses 
is further decreased to 1.23 eV, with an error or 0.13, 
incompatible with zero at high significance. 

Figure [4] shows the one dimensional marginalized pos- 
teriors of the joint analyses for TO4 in all the cases re- 
ported in Table |IIH together with the SBL posterior of 
the (3+1) analysis. Concerning the SBL posterior we can 
appreciate that the zero mass region is highly excluded 
and the maximum probability is characterized by three 
peaks between 1 eV and 2 eV, plus a small peak close to 
3 eV. In the combined analysis (which includes the SBL 
posterior as a prior), the joint posterior closely follows the 
SBL posterior, with basically the same best-fit located at 
7M4 = 1.27 eV. The only relevant difference of the joint 
posterior with respect to the SBL posterior is the sup- 
pression in the joint posterior of the SBL high-mass peak 
close to 3 eV, that is in fact highly disfavoured by cosmo- 
logical data, especially when the matter power spectrum 
is included in the analysis. We notice also that the inclu- 
sion of SBL information alleviates the tension between 
CMB and matter power spectrum: in fact, CMB-only 
data exhibit a best fit value close to 2 eV and a detec- 
tion of the single mass eigenstate 7714 = (1.72 ± 0.65) eV, 
while the inclusion of SDSS data shifts the results to- 
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Parameters 


WMAP9+SPT 


WMAP9+SPT 
+SBL 


WMAP9±SPT±SDSS 


WMAP9+SPT+SDSS 
+SBL 




0.02256 ± 0.00037 


0.02249 ± 0.00035 


0.02233 ± 0.00034 


0.02230 ± 0.00034 




0.131 ±0.008 


0.131 ±0.007 


0.127 ±0.007 


0.128 ±0.007 


Os 


1.0412 ±0.0011 


1.0411 ± 0.0011 


1.0411 ±0.0011 


1.0412 ±0.0010 


T 


0.083 ±0.013 


0.083 ±0.013 


0.081 ±0.012 


0.080 ±0.012 


Us 


0.959 ±0.011 


0.962 ± 0.009 


0.963 ±0.011 


0.958 ± 0.009 


log(10i°^.) 


3.222 ±0.041 


3.218 ±0.036 


3.212 ±0.036 


3.226 ±0.034 




0.65 ±0.22 


0.69 ±0.21 


< 0.96 


< 0.83 


7714 [e V J 


i . ( z zt U.DO 


1 97 -1-0 10 


<^ z.uy 


l.Zo zt U. lo 


Jfo[km/s/Mpc] 


68.6 ±2.1 


68.9 ± 1.7 


69.3 ± 1.9 


68.1 ± 1.4 




0.668 ±0.061 


0.692 ± 0.036 


0.766 ± 0.036 


0.744 ±0.034 




0.327 ±0.035 


0.325 ±0.030 


0.311 ±0.025 


0.324 ±0.027 


Xmin 


8274.1 


8274.3 


8326.4 


8327.5 



TABLE III. (3±1) analysis - Values of the cosmological parameters and their 68% confidence level intervals in the case of one 
additional massive sterile neutrino, with mass m4 and with multiplicity 7V4. The (3±1) SBL is applied where specified. 
Upper bounds are quoted at 95% C.L. 




FIG. 5. (3±1) analysis - Two dimensional marginalized 68% and 95% confidence level regions in the plane N4 - 7714 for the 
different combinations of datasets reported in Tables [Till 



wards lower values of 7714 and is consistent with a zero 
value of the mass eigenstate. 

As discussed above, the value of A^4 is almost uncon- 
strained when SDSS data are included, while in the case 
of CMB-only there is evidence for an extra massive ster- 
ile neutrino. This is manifest also in Figure [5] where the 
the two dimensional marginalized 68% and 95% C.L. re- 
gions in the plane (iV4,?Ti4) are plotted for the different 
combinations of datasets reported in Tables We can 
clearly see the effect of the SBL data: the constraints 
on the mass are strongly tightened but there is almost 
no effect on N4. Concerning the degeneracy between the 
number of sterile states and their mass, a negative corre- 



lation emerges when the matter power spectrum is taken 
into account and the result is a divergence of TO4 as 
is approaching zero. 



B. Bi— dimensional (3+2) Joint Analysis 

The results of the joint analysis for the (3+2) scheme 
are reported in Table |IV| and in Figure |6] and [Tj The 
analysis refers to the case of Eq. ( 12 1. 



Concerning the masses, the value of the heaviest mass 
eigenstate {m^) is always significantly deviating from 
zero: when CMB-only data are considered the effect is 
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Parameters 


WMAP9+SPT 


WMAP9+SPT 
+SBL 


WMAP9±SPT+SDSS 


WMAP9+SPT+SDSS 
+SBL 




0.02270 ± 0.00038 


0.02267 ± 0.00037 


0.02242 ± 0.00036 


0.02235 ± 0.00035 




0.139 ±0.010 


0.139 ±0.009 


0.133 ±0.008 


0.131 ±0.008 


Os 


1.0408 ±0.0011 


1.0409 ± 0.0011 


1.0406 ±0.0011 


1.0409 ±0.0011 


T 


0.083 ±0.013 


0.083 ±0.013 


0.082 ±0.013 


0.080 ±0.012 


Us 


0.961 ±0.012 


0.960 ±0.010 


0.967 ±0.012 


0.960 ±0.010 


log(10i°^.) 


3.228 ±0.041 


3.232 ±0.038 


3.213 ±0.037 


3.225 ±0.034 


N4 


< 0.92 


< 0.92 


< 0.92 


< 0.85 


1714 [e V] 


< 2.11 


1.20 ±0.30 


< 1.34 


0.95 ±0.30 




< 0.93 


< 0.91 


< 0.78 


< 0.62 


m5[evj 




i.yo =h U.4o 


i.40 ± U.D/ 


1 c^o 4- n /10 


J/o[km/s/Mpc] 


68.5 ±2.3 


67.9 ± 1.9 


69.6 ±2.1 


68.2 ± 1.5 




0.642 ± 0.060 


0.648 ±0.037 


0.752 ±0.036 


0.734 ±0.032 




0.346 ± 0.038 


0.352 ±0.037 


0.321 ±0.027 


0.331 ±0.026 


Xmin 


8274.1 


8274.8 


8326.4 


8327.8 



TABLE IV. (3±2) analysis - Values of the cosmological parameters and their 68% confidence level intervals in the case of two 
additional massive sterile neutrinos, with masses m4 and ms multiplicities A^4 and N5. The bi-dimensional (3+2) SBL is 
applied where specified. Upper bounds are quoted at 95% C.L. 




FIG. 6. (3±2) analysis - One dimensional marginalized posterior probabilities for 7714 (left panel) and (right panel). The 
thick (green) and thin (blue) lines refer to the case of Table IV with and without SBL prior, respectively. Solid lines stand 
for the analysis on CMB-only (WMAP9±SPT) data; dot-dashed lines refer to the inclusion of information coming from the 
matter power spectrum. The (red) dashed line shows the bi-dimensional 3±2 SBL posterior marginalized over m4 or ms. 95% 
C.L. upper bounds on the masses are reported as vertical lines. 



at the 3.5cr when the SBL are not included, and grovirs 
to a 4(7 effect when SBL priors are considered. In both 
cases, the neutrino mass is around 2 eV. When large- 
scale-structure data are included, the preferred value for 
the mass decreases around 1.5 eV, with a significance of 
2.1(3.2)cr without (with) the inclusion of SBL data. SBL 
results, which have a clear preference for massive neutri- 
nos close to 1.6 eV, reinforce the cosmological results and 



therefore induce a clear increase in the confidence for the 
mass determination. This fact occurs also for the lighter 
mass eigenstate (m^). While cosmological data do not 
require the lightest neutrino to be massive (bounds are 
at 2.11 eV for CMB-only data and 1.34 eV when SDSS 
is included), SBL priors induce a clear preference also for 
non-zero lightest mass eigenstate around 1 eV (an effect 
close to the 4tT level). 
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FIG. 7. (3+2) analysis -Two dimensional marginalized 68% and 95% confidence level regions in the plane (m4 — m^) for the 
different combinations of datasets reported in Table [rv] The two dimensional marginalized 95% SBL posterior probability of 
the 3+2 model is also overplotted (red solid line). 



The multiplicity parameter of each mass eigenstate 
is basically unconstrained for all the different combi- 
nations of datasets. Full and standard thermalization 
{Ni = 1,1 — 3,4) is not allowed, but fractional occu- 
pations as large as 0.92 are possible. A notable excep- 
tion occurs for the heaviest eigenstate {m^) when SDSS 
data are included: in this case, the multiplicity param- 
eter needs to be smaller, not exceeding 0.78 for the full 
combination of cosmological data and further decreasing 
to 0.62 when the SBL information is included. From the 



last column of Table |IV] we can conclude that cosmolog- 
ical and SBL data are consistent with the existence of 2 
extra sterile neutrinos, provided that they are not fully 
thermalised. 

The same results can be appreciated with more level 
of detail by directly looking at the behaviour of poste- 
rior probabilities. Figure |6] shows the one dimensional 
marginalized posterior for (left panel) and (right 
panel) in all the cases reported in Table IV together 
with the SBL posteriors of the (3+2) scheme, marginal- 
ized over the parameters not shown in the figure. The 
SBL posteriors are characterized by the clear exclusion of 
zero values for both and ms . The cosmological poste- 
riors of TO4 prefer a zero mass value, even more so when 
the matter power spectrum information is included in the 
analysis. When the SBL information is applied, the 
null values for the neutrino masses become clearly dis- 
favoured. The matter power spectrum effect turns out in 
a shift of the preferred mass range towards lower values. 

Figure [t] shows the two dimensional marginalized 68% 
and 95% confidence regions in the plane 7714 - for 
the different combinations of datasets reported in Table 



IV plus the two dimensional marginalized 95% SBL pos- 
terior of the 3+2 model. Even when the matter power 
spectrum is included in the cosmological analysis, the 
blue contours referring to cosmological-data only show 
a good agreement with the red 2a confidence region ob- 
tained by the SBL analysis. 



C. Reduced (3+2) Joint Analysis 

Finally, in Table |V] and in Figure [8] and [9] we show the 
results of the join (cosmological + SBL) analysis in the 
(3+2) scheme, where we have post-processed the chains 
of the MCMC to get constraints only on the sum of the 
mass eigenstates and on the total number of extra species, 
by applying a SBL prior on (777.4 + 7775) obtained by real- 
izing a one-dimensional SBL posterior on the sum of the 
two masses. 

Interestingly in this case the sum of the neutrino 
masses always shows a clear preference for a non zero 
value, even when only cosmological data are considered: 
the confidence level for CMB-only data is 3.1(7, and de- 
grades to 2.2(7 when SDSS data are included. The inclu- 
sion of the SBL information sizably strengthen the result, 
driving the confidence of non-zero masses close to a 4ct 
level. 

Figure [8] shows these behaviours by reporting the one 
dimensional marginalized posterior for 7714+7775 in all the 
cases reported in Table |Vj together with the posteriors 
obtained by analyzing the SBL data within the frame- 
work of the (3+2) model and marginalizing over 7774+7715. 
The situation is analogous to the one we found in Section 
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Parameters 


WMAP9+SPT 


WMAP9+SPT 
+SBL 


WMAP9±SPT±SDSS 


WMAP9+SPT+SDSS 
+SBL 




0.02270 ± 0.00038 


0.02270 ± 0.00039 


0.02242 ± 0.00036 


0.02241 ±0.00035 




0.139 ±0.010 


0.139 ±0.009 


0.133 ±0.008 


0.132 ±0.008 


Os 


1.0408 ±0.0011 


1.0407 ± 0.0011 


1.0406 ±0.0011 


1.0407 ±0.0011 


T 


0.083 ±0.013 


0.084 ±0.013 


0.082 ±0.013 


0.081 ±0.012 


Us 


0.961 ±0.012 


0.963 ±0.011 


0.966 ±0.012 


0.964 ±0.010 


log(10i°^.) 


3.228 ±0.041 


3.225 ± 0.040 


3.213 ±0.037 


3.217 ±0.035 




0.97 ±0.31 


1.00 ±0.31 


0.75 ±0.32 


0.69 ±0.29 




Q f\A -\- n 07 


9 Sin -1- n 71 


i.yo zt u.oy 


z.oo zt U.Dl 


i7o[km/s/Mpc] 


68.5 ±2.3 


68.7 ± 2.2 


69.6 ±2.1 


69.1 ± 1.8 




0.642 ± 0.060 


0.647 ±0.054 


0.752 ±0.036 


0.744 ±0.033 




0.346 ± 0.038 


0.344 ±0.037 


0.321 ±0.027 


0.325 ±0.026 


Xmin 


8274.1 


8275.1 


8326.4 


8327.6 



TABLE V. Reduced (3±2) analysis - Values of the cosmological parameters and their 68% confidence level intervals, for the 
data analyzed in terms of of the sum of the sterile neutrino mass eigenstates (m4 ± ms) and of their total effective multiplivity 
A^4 + A^s. The one-dimensional (3±2) SBL is applied where specified. Upper bounds are quoted at 95% C.L. 



1.0 



0.8 



> 0.6 



sbl 

cosmology 
cosmology+sbl 




(m^+m^) [eV] 



FIG. 8. Reduced (3±2) analysis - One dimensional marginal- 
ized posterior probabilities for (m4 ± ms). The thick (green) 
and thin (blue) lines refer to the case of Table |V] with and 
without SBL prior, respectively. Solid lines stand for the anal- 
ysis of CMB-only (WMAP9±SPT) data; dot-dashed refer to 
the inclusion of information from the matter power spectrum. 
The (red) dashed line shows the one-dimensional 3±2 SBL 
posterior on (m4 ± ms). 95% C.L. upper bounds are also 
reported (vertical lines). 



|VB[ The cosmological posteriors alone show a preference 
for a non-zero value for the sum of the masses. As usual, 
when SDSS data are considered, the posterior is shifted 
towards lower values for the sum of the masses, but the 
preference remains for masses around 1 eV. Concerning 
the SBL posterior, the maximum probability on the sum 



of the masses is around 2 eV as expected, with a tail at 
higher masses. This higher mass tail is suppressed in the 
joint posterior, when the cosmological analysis considers 
also the matter power spectrum information. Instead, a 
higher mass tail can be seen in the joint posterior when 
the combined analysis takes into account CMB-only data 
as cosmological datasets. 

Figure [9] shows the two dimensional marginalized 68% 
and 95% confidence level regions in the plane {N4 + A^s) 
vs. (m4 ± for the different combinations of datasets 
reported in Tables |Vj Once again, we can clearly see the 
effect of the SBL data: the constraints on the sum of 
the masses are tightened but there is almost no effect 
on the number of extra massive sterile neutrino species. 
Furthermore, in Figure |9] we recognize the same degen- 
eracy between the sum of the mass eigenstates and their 
multiplicity that was seen in Figure [5] As we have al- 
ready discussed, this degeneracy appears when matter 
power spectrum information is taken into account. In- 
terestingly, here this degeneracy remains also when the 
SBL posterior is applied. 



VI. CONCLUSIONS 

The case for Extra Dark Radiation in cosmology has 
recently become very complex and somewhat controver- 
sial: the new WMAP9 and SPT data are still pointing 
toward an extra component while the new ACT results 
became compatible with the standard cosmological value 
Ncs — 3.046. On the other hand short-baseline neu- 
trino oscillation experiments keep confirming the LSND 
anomaly and the presence of one or two sterile neutrinos. 

In this framework, our analysis provides an update of 
the cosmological results in the sterile neutrino context, 
taking into account all the present CMB data (WMAP9, 
SPT, ACT) and investigating all the possible parame- 
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terizations (massless; massive; 3 active massless plus a 
varying number of massive sterile states). 

We find that even in the context of such very gen- 
eral models the ACT and SPT data are not compatible, 
and that the results are strongly affected by this discrep- 
ancy. This in turn affects the inferred neutrino mixing 
and thermalization parameters. The discrepancy is to 
some extent alleviated when BAO data plus a prior on 
Hq from the HST analysis are included: in this case the 
effective number of (massive) neutrinos is 3.60±0.35 with 
a total effective mass 0.59 eV. 

Given this discrepancy we have not used the ACT 
dataset in the rest of our analysis (because it also leads 
to spuriously stringent bounds on the neutrino mass), 
and we considered only WMAP9 and SPT in combina- 
tion with SDSS-DR7. In order to analyze the cosmo- 
logical evidences in the context of neutrino mass mod- 
els and short-baseline neutrino oscillation data, we have 
used the results from the SBL as a prior in the analyses 
of cosmological data. We specifically considered models 
with one or two extra sterile neutrinos (denoted by (3-1-1) 
and (3-1-2) schemes): we first performed a full analysis of 
the SBL data and then we used the SBL posterior prob- 
abilities on the sterile neutrino masses as priors in the 
MCMC analysis of the cosmological datasets. 

We found that the inclusion of the SBL priors induce 
tight bounds on the sterile neutrino mass eigenstates, 
and mildly constrain the fractional contribution of the 
extra neutrinos to the total energy budget. In the (3-1-1) 
scheme we obtain m4 = (1.27 ± 0.12) eV when CMB- 
only data are considered as cosmological datasets, and 
m4 = (1.23±0.13) eV when SDSS-DR7 data are also in- 
cluded. We also notice that there is an evidence for a non 



zero value of the single mass eigenstate (although with a 
larger uncertainty) even without SBL priors when only 
CMB data are considered: m4 = (1.72 ±0.65) eV. In the 
(3-1-1) scheme, the inclusion of SBL information does not 
significantly constrain the multiplicity parameter, which 
could be as large as 0.96. Instead, if SDSS data are not 
included in the analysis, the multiplicity parameter can 
deviate from zero at about 3(t. 



In the (3-1-2) context, the inclusion of SBL information 
induces relatively tight intervals for the two mass eigen- 
states: mi = (1.20 ± 0.30) eV and = (1.96 ± 0.48) 
eV for CMB-only data; m4 = (0.95 ± 0.30) eV and 
TO5 = (1.59 ±0.49) eV when SDSS-DR7 is also included. 
In this scheme, the multiplicity parameters are again only 
slightly bounded from above, except for the heaviest mass 
eigenstate: when SDSS data are included, cannot ex- 
ceed 0.62, meaning that the sterile neutrino need to be 
only partially contributing to the energy density of the 
Universe. 



In conclusion, we found that the SBL data exhibit 
a good agreement with the new updated cosmologi- 
cal bounds on neutrino masses: this occurs both when 
considering CMB-only data (WMAP+SPT) and when 
adding information from SDSS. Either one or two extra 
sterile neutrinos, with a mass at (or even above) the eV 
scale (as dictated by the SBL results) and not fully con- 
tributing to the energy density (as can occur, e.g., in the 
case of partial thermalization) are therefore fully com- 
patible with current cosmological measurements. 
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